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Abstract 

A relativistic quark model is extended to incorporate chiral and gauge symme- 
tries. We obtain the DD*Tr and DD*^ couplings and find that the ratio T{D*^ — > 
D^TT^)/T{D*^ — > D^^) constrains the charm-quark mass close to 1.45 GeV. Large 
l/rUc corrections appear in the heavy-quark contribution to the DD*'y coupling. The 
model is extended further to describe seven excited D meson states. We find that 
semileptonic B decays into the ground and excited D meson states do not account for 
the total semileptonic decay width of the B^. The nonresonant contributions to the 
processes appear to be large enough to account for the discrepancy. 
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Introduction 



We have recently developed a fully relativistic model of hadrons containing one heavy quark 
|I|]. Hadronic matrix elements are represented by quark loop graphs with momentum depen- 
dent interaction vertices between hadron, heavy quark and light quark. The only parameters 
in the model are the quark masses. We previously applied the model to B ^ D^*Ht7 ||, 
B K*-f §, Afe AJT^ §, and Q'C^ Q^;Ut7 §. The model does not rely on any expan- 
sion in inverse powers of heavy-quark masses mq. However, QCD imposes severe constraints 
on the form of the heavy-quark expansion, and not all models automatically satisfy these 
constraints 0. We have verified that the heavy-quark expansion of our model does satisfy 
these constraints. 

Confinement is modeled in a simple way; we simply drop the imaginary parts of our loop 
graphs. This prescription preserves the symmetry structure of the theory, and our hypothesis 
is that it may give reasonable results as long as external momenta of a diagram are such as to 
be far removed from whatever unphysical singularites the diagram possesses. (A discussion 
of the singularity structure of our three-point functions, including anomalous thresholds, 
may be found in [^.) We are finding that this simple, most naive model of confinement 
is yielding some very encouraging results. It also acts as a useful preview to the general 
Bethe-Salpeter approach. 

Although our results are based only on a model, they suggest that some caution is required 
in the use of more systematic attempts to extract information from QCD. In particular, some 
important results in the literature depend on an expansion in powers of 1/ rric. We question 
whether this expansion is sufficiently understood and under control at present. Our model 
provides an explicit example of how far from the heavy quark limit the charm mass could, 
in fact, be. We will later discuss an example of a quantity receiving large l/rric corrections, 
and which also contributes to a measured branching ratio. 

One of the prime applications of the heavy quark expansion is in the extraction of \Vcb\ 
from the processes B — > D^*Hv and A^ — > Kctu. In the heavy-quark limit oo, 
these processes are normalized in a model-independent way at zero hadronic recoil. But for 
any realistic in our model there are always large positive deviations in at least one of 
B — s> D*£v 01 B ^ Dtp p|, concurrently with rather small deviations in A5 h^ctP 0. 
The large deviations are a sensitive function of the quark masses, and their origin lies in the 
hyperfine interaction which splits the ground-state pseudoscalar and vector D mesons. This 
would lead to a discrepancy between the values of extracted from these various processes 
under the assumption that the deviations from the heavy-quark limit are negligible. 

In Fig. 1 we show that our model's prediction for the shape of the decay spectrum of 
B D*iT' agrees well with the most recent CLEO data 0. If we ignored our correction to 
the overall normalization and forced the heavy-quark limit normalization at zero recoil we 
would obtain |V"cb| = 0.036. Our correction to the normalization implies that the actual \ Vcb\ 
could be 10% or 15% smaller. 

This paper is divided into three parts. In Parts I and II we describe two major extensions 
of the model. In Part I we incorporate chiral and gauge symmetries to deal with various 
processes involving pions and photons. In particular the ratio r{D*^ D^7T^)/r{D*^ 
D^'-y) is well measured experimentally and provides a sensitive probe to the physics of the 
light degrees of freedom in the D mesons. This will help constrain the only parameters 
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appearing in our model, the quark masses. Our model reproduces this branching ratio for 
a choice of quark masses very close to what we have previously used, and in particular it 
contrains our model rric to be close to 1.45 GeV. We have observed before that rrib — uic was 
much more constrained by our model than rrib+mc] the constraint reads mb—rric = 3.36±0.03 
GeV (see the wedge-shaped region in Fig. 6 of 0). With our new constraint on nic, our 
model ruh must be close to our previously-used value of 4.8 GeV. 

One ingredient in our calculation of r(L)*° — »• D'^n^)/r{D*^ —>■ D^'^) is the heavy-quark 
contribution to the D*D'~f coupling. We find large l/rric corrections to this quantity, which 
are analogous to the large corrections we find in semileptonic B decays. For the D*Dtt and 
D*D'y couplings we find values which are one-third as large as the nonrelativistic quark model 
predictions. The small D*D7i coupling has implications for calculations in heavy-quark chiral 
perturbation theory, and in particular the quantity Ad = {Mjj* — M/jJ — (M^* — Mr,). We 
also include a discussion of charge radii. 

We then consider the nonresonant contributions to the processes B D(*)7r£z7 in the 
context of the problem of the missing exclusive semileptonic B decay modes. The branching 
fraction for B £z7+ anything is (9.5±1.6)%, while those for the dominant exclusive modes 
B° D+iVandB° D*(2010)+£z7 are (1.9±0.5)% and (4.4±0.4)%, respectively [0. The 
model indicates that the nonresonant contributions play an important role in accounting for 
the discrepancy. 

In Part II we extend the model to include seven excited charmed mesons. From the 
masses (or estimates of the masses) of these mesons we are able to determine appropriate 
model vertex functions. We calculate and display the various amplitudes for semileptonic 
B meson decays into these excited states. Our results for the various decay widths and 
electron energy spectra are in qualitative agreement with those of other authors. We confirm 
in particular that the decay rates to excited Z)'s are too small to account for the above 
mentioned discrepancy in semileptonic B decays. 

In Part III we present some details of the heavy-quark limit of the model as it applies to 
the excited D's. The full-model form factors deviate considerably from their values in the 
heavy-quark limit, as expected since the masses of these states are significantly greater than 
the charm-quark mass. We include the heavy-quark-limit results primarily for comparison 
with other models. 

Also in Part III we study two modifications of the model. We first study the possible 
modification necessary to describe strange mesons, and in particular the process B^ 
K e'^v. We address the question of whether the s should be treated as heavy or light. 
And second, we investigate the effect of modifying the form of the momentum damping at 
the basic quark-antiquark-meson vertex of the model. We consider an exponential damping 
factor, and we recover our previous conclusion that there is always a large positive deviation 
from the heavy-quark limit in one oi B ^ D*iT7 or B DiV. 
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Part I 



Gauge and chiral symmetries 



We begin with the hght-quark triplet and triplets of pseudoscalar or vector meson fields Hq 
with heavy quark Q. 



Q Hq ^ {Q%Qd,Qs) 



{B-,B\b''^) and q 



iD'',D+,Dt) 



( 

d 



Under local SU(3)yxyi transformations with external vector and axial-vector fields 
and Afj^ transforming as 



(2) 



the quark fields transform according to Q ^ Q and q — > Gq, where G = exp{i[9v + OAjb])- 
We define ^ = eyop{iM^5/ fn), where ~ 132 MeV and 



M 



K+ \ 



V 



TT 

K- 



K 



75^ / 



Then ^ transforms as ^ — > gr^G^ 

The gauge-invariant version of the interaction Lagrangian is 



Ant = Y^Y^j dy Q{x)CH{x,y)FH{x - y) + h.c., 



Q=b,c H 

where the function Fh{x — y) is the Fourier transform of the damping factor, 



For the ground-state mesons, n — 1 and Ch{x, y) are given by 

H Cnix^y) 



1 0^/2 -iHQ{xYi^K{x,y)(,{y)q{y) 
1 1-/2 iflQ{:x)K{x,y)^{y)q{y). 



The quantity K{x, y) is a path-ordered exponential 

K{x,y) = Vex^i^if^dz^T^{z)Y 

where 



(3) 



(4) 



(5) 



(6) 



(7) 
(8) 



Under the SU(3)v'xA transformations, K{x,y) g{x)K{x,y)g\y), and the meson fields 
transform as Hq Hqq'^ . 

Mass terms involve an SU(3)vxA-invariant piece and an SU(3)yxA-breaking piece. 

>^^mass = + Mq)q, (9) 

where we take m ~ 250 MeV and Mq = diag{ml^" , mf', mf"'') ~ (0,0, 170) MeV. When 
A^q = 0, the full Lagrangian 

C = Q{i^ - mQ)Q + q{i^ + f+ 475)g + £mass + Ci^t (10) 

is invariant under the local SU(3)yxA transformations. 

D* Dtt and D* Dj 

In this section we describe our calculation of r(D*° D^7i^)/r{D*^ — * -0%). This quantity 
depends sensitively on the only two parameters appearing in our model for D mesons, the 
charm mass and the light quark mass m. We will show that a parameter choice which gives 
an acceptable result is rric = 1450 MeV and m = 270 MeV. The main point is that these 
values are close to the values, rric = 1440 MeV and m = 250 MeV, appearing in our previous 
work. For consistency with our previous results we will continue to use our old parameter 
choice in all other sections of this paper. Except where noted, our results in other sections 
do not display such a large sensitivity to quark masses. 

In this section we shall also be comparing the D meson results to B meson and heavy- 
quark limit results. For the b quark mass we adopt our previous value, rUb = 4.8 GeV, as 
explained above. 

We first consider the matrix element 

U{H%M'v')n+{l)\H*+{Mv,e)) = ^MM^{u + l)h^^{u)l ■ e. (11) 

We define Qu'R-k = h^^{l). This can be calculated using the pion couplings derived above, 
or it can be extracted more simply by inserting the axial current on the light-quark line. 
Our results are gD*Dn = 0.28 and gB*B-K = 0.32. The first two terms in the heavy-quark 
expansion are 

gH*Hn = 0.3A-0.17—, (12) 

where A = M — itlq ~ 500 MeV is the lowest-order difference between the meson and heavy- 
quark masses. These numbers are in agreement with other relativistic models and sum rules 
|T^|T^[|I^, and they are in marked contrast to the nonrelativistic quark model which gives 
a value of unity. 

We next consider the matrix element of the electromagnetic current, and we distinguish 
the light- and heavy-quark contributions. 

{HQ{M'v')\j;-^-\H*Q{Mv,e)) = VMArie^hf^u;) + eQh^^uj)]te,.,^e%'Pv- . (13) 
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For the hght-quark contribution we find h^'^{l) = 1.56 and h^"'{l) = 5.00. The first two 
terms in the heavy-quark expansion are 

/if^(l) = 0.50^ (^1 + 0.47^^ . (14) 

The hght-quark contribution is not predicted by heavy-quark symmetry, but the nonrela- 
tivistic quark model predicts a value for h^^{l) of order the heavy-quark to light-quark mass 
ratio. Our result is about a factor of three smaller, which mirrors the situation with the pion 
coupling. In both cases, the relevant diagram involves either a vector or axial coupling to the 
light quark line. The suppression associated with the relativistic nature of the diagram seems 
to apply to both cases equally. We also note that in both cases the heavy-quark expansion 
provides a good representation of the full-model results. 

The situation is different for the heavy-quark contribution to the DD*'~f couplings. We 
find hg'^ll) = 1.54 and /;.q'^(1) = 1.16, to be compared with the first two terms in the 
heavy-quark expansion 

^?"(1) = 1 + — (1-6(1)). (15) 
niQ 

The latter gives hQ''{l) = 1.35 and /iQ^(l) = 1.10. (^3(1) is negligible in our model [|l|.) The 
deviations from the heavy-quark expansion are much more significant here than they are for 
the pion and photon couplings to the light quark. As we have said, these large corrections 
are analogous to the large l/rric corrections we find for Vet- 

We now use the above results to extract the Dtt and Dj decay widths of the D*. We find 
T{D*+ D°7r+) + T{D*+ D+Ti^) = 20.9 KeV and T{D*'^ D^tt'^) = 9.4 KeV. For the 
radiative decay branching fractions we also display the effect of the leading SU(3) breaking 
corrections calculated in |Tl|. We obtain 



experiment with corrections without corrections 
B (D*o ^ DO7) (%) 36.4 ±2.8 37.8 45.6 ^^^^ 

B ^ D+7) (%) 1.1 7 2.9 2.1 

The agreement in the first two columns is a refiection of the quark masses we have chosen. 

We consider briefiy the implications for the S'?7(3)-violating quantities Ah = {Mh* — 
MhJ — {Mh* — Mh). It was stressed in []T3[ that the measured values, A/j = 0.9 ± 1.9 MeV 



and As = 1.2 ± 2.7 MeV, are much smaller than expected from the leading SU(3)-breaking 
one-loop contribution. It was estimated that this source alone gives A^'^ ^ —47 MeV 
and /^^'^ '°°P ~ _X5 assuming that g]:,,£,T, = 0.5. Using our values of go'Dw and gB*B-K 
instead would yield the less problematic numbers A^*^ pa —7 MeV and A^° '""'^ pa —3. 

Charge radii 

Gauge-invariance of the Lagrangian (|TD|) implies extra Feynman rules for the matrix elements 
of the light-quark vector current between two mesons. The photon can attach to the nonlocal 
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vertices as well as to the quark lines. The new Feynman rule for a vertex with outgoing 
photon momentum q and outgoing light-quark momentum k is 



[-F + A2][-(A; + g)2+A2]' 

This must be included in the computation of the light-quark contribution to the electromag- 
netic charge radius of the meson (although it did not appear in the above calculation of the 
-D-D*7 coupling). 

Define the matrix element of the electromagnetic current by 

(^Q(p')l^ri^Q(P)) = [-^M<f) + <^QFQ{q')]{p + p'),, (18) 

where Cc = e„ = 2/3 and 6^ = 6^ = —1/3 and we have decomposed the form factor 
into light-quark and heavy-quark contributions with -Fq(O) = 1 = Fq(0). We define rj = 
[69Fj/(9g^(0)]^/^ for i = q,Q. As an example we compare the contributions to the charge 
radius of the ground state pseudoscalar {L = 0) D to one of the excited states treated in 
Part II, the scalar (L = 1) D^. 

H l/rg(MeV) l/rQ(MeV) 

D 300 830 (19) 

D* 690 830 

The 300 MeV vs. 690 MeV shows that the excited state is more compact than the ground 
state, as might be expected due to the higher energy present in the light-quark system in 
the excited state. This effect is not present in models of nonrelativistic quarks. 

We also note that the large light-quark charge radius of the D reflects a rapid damping of 
the form factor in q^ when the current attaches to the light-quark line. We will use this fact 
when discussing nonresonant contributions to semileptonic decays in which a pion attaches 
to the light-quark line. 



Nonresonant processes 

The data indicate that the decays 5^ D^tu and 5^ D*(2010)+£l7 contribute only 
66% of the total semileptonic width B ^ tP + anything. We will find in Part II that 
resonant decays to the eight lowest lying excited D's account for a further 6%. This leaves 
28% missing. It is safe to assume that this will not be made up by decays into even higher 
lying D excitations. And neither will it be decays to final states other than charm, in view 
of the smallness of \Vub\- The implication is that a significant portion of the total width 
must be into nonresonant charmed final states. The nonresonant width must be of order 
1.6 X W^\Vcb\^ s-^ 

We consider nonresonant contributions to the processes 5^ D^*^7iiV in our model. The 
Lagrangian ([T0|) gives rise to three graphs, one with a pion coming off the light-quark line 
(proportional to the common mass m) and two with pions coming off the vertices. The sum 
of these three graphs is the same as an insertion of /~^/^7'^75 on the light-quark line, where 
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/ is the pion momentum. Because of the difficulties involved in computing the resulting 
four-point graph, we were able only to obtain the hadronic matrix element in the soft-pion 
limit, i.e. keeping only terms linear in I. In this limit, we define form factors qi{uj) and ri{uj) 
(with u = V ■ v' as usual) by 

UD%MDv')n+il)\c^,,{l - ^,)h\B\MBv)) = (20) 
qvi£f,upav''Fv"' - qAjf, - [g^a^ ■ v + qaJ ■ v']v^ - [qaJ ■ v + qaJ ■ v']v' . 



and 



MD*%MD-^v')7c+{l)\c^^{l-^5)b\B\MBv)) = 

e ■ v[rv^l^ + {ry^v ■ I + ty^v' ■ l)v^ + {ry^v ■ I + ry^v' ■ l)v'^ 
+ {rveV ■ I + ry^v' ■ l)e^ + {ry^v^ + ry^v'^)l ■ e 

-le^up.eyA.v'Pv'' + I'irA.V' + rA.v"")]. (21) 



The values of the above form factors at = 1 are 

qy QAi qA2 QA3 qA4 QAs 



14.8 24.0 -3.5 -8.2 -4.6 -2.7 

and 

ry, ry^ + ry^ + ry^ + ry^ ry^ + ry, ry^ + ry^ ta^ + ta^ 
-1.7 11 -0.1 -20.1 19.5. 



(22) 



(23) 



The rather large magnitude of some of these dimensionless hadronic quantities is striking. 
The amplitude vanishes for vanishing pion momenta, and it rapidly increases for increasing 
pion momentum up to some momenta of order 300 or 400 MeV (roughly the inverse of the 
charge radius), after which the amplitude will decrease. We choose a small value of the pion 
momentum (equal to the pion mass) well below where the amplitude peaks. We set v = v', 
let the massless electron and antineutrino emerge with equal energies, and then perform the 
four-body phase space integral with this constant matrix element squared. The result is 
T(B^ D^TT+eu) ~ 0.6 X 10^^\Vcb\^ s-^ and D*^TT+iV) ~ 0.3 x 10^^\Vcb\^ s'^Q The 

decays with charged D and neutral pion are half as large. The total nonresonant width is 
then ~ 1.8 x 10^^|Vcfcp s~^, in rough agreement with the value deduced above. While this 
crude estimate is far from being a definite prediction, it does make reasonable the possibility 
that nonresonant decays make an important contribution to the semileptonic width of the 
B meson. 



Part II 

Description of model for excited mesons 

Our notation is J^^, where labels the radial excitation, is the meson spin-parity, and 
ji is the total angular momentum of the light degrees of freedom. It is the sum of orbital 

-'^ These results are also quite sensitive to the quark masses used. For example, if our "new" values of the 
masses are used then T{B D'^n^tU) increases by 40%. 
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angular momentum L and light-quark spin si = 1/2. This paper will treat the eight 



states shown in (p^ , where the masses in the last column corresponding to question marks 



are estimates. The naming convention is that of the Particle Data Group [IC] 

L ji doublet name mass 



1 1 
1 2 




D 1869 

D*(2010) 2010 

D*{7) 2465 

D^{7) 2270 (24) 

Di(2420) 2421 

D*l2AQ0) 2465 

D*{1) 2800 

Dsf? 2800 



In addition, we will consider the first radial excitation 2 Oj^y, with an estimated mass of 2440 



MeV. As is commonly done |T^, we will ignore possible mixing between 1 and 1 Igy,- 

The interaction vertices between hadron, heavy quark, and light degrees of freedom are 
chosen in the following way. Each vertex is the product of some gamma-matrix structure 
and a damping factor. The gamma-matrix structure is chosen to be identically equal to the 
form determined by heavy-quark symmetry [|l^]. This is required for consistency with QCD 



in the heavy-quark limit. The damping factor is chosen to be a sum of terms of the form 
[Z'^/ (— A;^ + A^)]", where the least power n is chosen to be the smallest which guarantees the 
convergence of the relevant integrals, and where Z and A are different for each state. We 
use standard quark propagators throughout. 

We take the interaction Lagrangian to be that given in where Ch{x, y) and the powers 
n are given by 



H 






n 




-iH{x)-f5q{y) 




1 


10+2 

11^/2 


i^{x)q{y) 
iH{x)q{y) 

ii^{xh5q{y) 




1 
1 
1 


1 I3/2 

1 9+ 

^3/2 


lm-{x)^^d,q{y) 


')]d.q{y) 


2 
2 


1 I3/2 

1 23/2 

2 0^/2 


lm^{x)^,^,d,q{y) 
-iH{x)-i^q{y) 


)}duq{y) 


2 
2 

1 + 2. 



(25) 



Wherever a velocity v occurs in (|2^), it is the velocity of the external meson (and not that 
of the heavy quark) . 

The only parameters of the model are the quark masses; the various Z's and A's are fixed 
via the meson "mass functions" S(p^). These are defined in terms of the meson self-energy 
graphs, which are given by zS(p^), —ig^yT^ijP') + ■ ■ ■ and ig^j.^gpaT^ijP') + ■ ■ ■ when J = 0, 1 
and 2, respectively. First A is fixed by setting S(M^) = 0. Then, Z is fixed by setting 
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S'(M2) = 1, where prime denotes differentiation with respect to p^. This latter condition is 
equivalent to conservation of the heavy-quark vector current by the Ward identity. 

The situation is a little more complicated in the case of the radial excitation 2 0^^2 of ^^e 
ground-state pseudoscalar 1 O^^g- take the damping factor to be a linear combination of 
n = 1 and n = 2, which we write as 



-P + A2 



1 + / 



A2 



-P + A2 



(26) 



First, / and Z are determined as functions of A at fixed quark and meson masses by the zero 
and the slope of the mass function. Then, A is determined by the vanishing of the matrix 
element (1 0]~y2lQ7At'5|2 0^/2) 'f — 0; required by current conservation. 

The values of the various constants are given in (|27|) in MeV. Quark masses rric = 1440 
and m„ = 250 MeV have been used. 



H 



A 



H 



1 

1 1 

1 

1 1 

1 1 

1 2 

1 1 

1 2 

2 



1/2 

1/2 

+ 

1/2 

+ 

1/2 
+ 
3/2 

+ 

3/2 
3/2 
3/2 
1/2 



537 

762 

1217 

1065 

1400 

1359 

1976 

1745 

1430 



745 

1193 

1838 

1881 

2143 

2122 

3231 

2804 

2594 



(27) 



We find / = —1.1 for the 2 0^/2 state. In addition, we need the values appropriate to the B 
and B*: mi, = 4800, Ab = 621, Ab- = 702, Zb = 1103, and Zb- = 1290 MeV. 



The resonant processes B Hip 

We now consider the hadronic matrix elements 

-Mf = -^j^{H{M'v')\ci,{l-imB{Mv)), (28) 
where H denotes any of the D states. The matrix elements may be written as 
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H 

llV2 



+ 



/2 



1 0^ 

1 1 + 

-'-3/2 
1 ^3/2 



1 1 
1 2 



3/2 



3/2 

2 07/2 



Ml 

-5(c+,c_)^ 

—T{dA, dvi, dv2, dv3)fiu£*'^ 

—T{fA, fvi, fv2J fvs)fiu£*'' 

T{gv, 9Ai , 9A2 , 9A,)t^ue*'"^vx 

T{hv, flA^, hA2, ^As)/^!^^*^ 

~T{kA, kvi, kv2, ky^) V\ 



(29) 



where 



S'(a+,a_)^ = a+{v + v')^ + a^{v - v')^ 



T{bv, bA^,bA2,bAs)ij,u 



ibv£^iupuv'''v'' - bA^9tiiy - (^As^^m + bA^v'^v^. 



(30) 



We take the form factors to be functions of u; = f ■ f e.g. a+ = a+(co'). 

The form factors a± and 6j of ( PP| ) and (|5D|) become proportional to the "original" Isgur- 

Wise function ^ = 'C1/2 iii the heavy-quark limit ^ 00. Similarly, c± and become 

proportional to ri/2, /« and gi to r3/2, and hi and fcj to ^3/2. The form factors l± for decay into 

(2) 

the radially- excited state become proportional to C\ i- The coefficients of proportionality 



are shown in (|31|). 



units 


form factors in given units 


6/2 


a+ = 1 
V = 1 
feA, = 


a„ = 

bAi = UJ + 1 

bA, = -1 


T"l/2 


c+ = 

c/vi = 2(cj - 1) 
dy^ = —2 


c_ = -2 
(ivj = 
dA = 2 


T"3/2 


f UJ^~\ 

/V2 - JVz - 

9A2 = 


f — 3 

fif^i = V^iu; + 1) 
5fA3 = -\/3 


6/2 


'^y -172 

ky^ = V5{UJ - 1) 
ky^ = —a/3 


hA2 + ^Ag = 

ky2 = 

^A = V3 


f(2) 
Sl/2 


= u; — 1 


L = 



(31) 



Our results for the various full-model form factors for B decays into the excited D states 
are represented in Figs. 2-4. In order to facilitate comparison with the heavy-quark limit, we 
have chosen not to plot the full-model form factors directly, but rather to divide out the co- 
efficients in (pTI). This puts all form factors on a comparable scale, and shows the "effective" 
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Isgur-Wise functions which would be extracted from the data under the assumption that 
the deviations from the heavy-quark hmit were small. In those cases where the coefficient 
vanishes at a; = 1, we have subtracted off the zero- recoil value before dividing by the coeffi- 
cient. For example, the curve labelled "ciy/' in Fig. 2 is equal to [dvi{uj) — dviil)]/^^^ — 1). 
In those cases where the coefficient vanishes identically, we have omitted the corresponding 
full-model form factor for brevity. 

Fig. 5 shows a comparison between the full-model form factors for B decays into the 
ground state pseudoscalar D and its first radial excitation. We note the strong suppression 
of the form factors for the latter. 

The various decay widths are given by 

th = iK^r 3 / ci^v/^T^^nM, (32) 

where r = M' /M. The functions J^h are given by 



IO1/2 Fo{a+,a_) 

1 1-/2 2FT(6y, bA,) + FLibA„bA,,bA,) 

10+2 ^o(c+,c_) 

1 I1/2 '^Ft^cIa, dvi) + Fi^dvi, dv2, dvs) 

1 I3/2 "^FxifA, fvi) + FLifvi, fv2^ fvs) 

1 23/2 (uj^ - l){Fr(c/y, c/Ai) + iFiigAi, 9A2, 9Ai)} 

1 I3/2 '^Fx^hv, hA^) + Fl^Haj^, hA2, hA^) 

1 2-/2 {uj^ - 1){Ft(A;a, hy,) + lFL{kv„kv2, ky,)} 

2 0^/2 Fo (/+,/_), 



where 



(33) 



Fo(a+,a_) = (^2-l)[(l + r)a+-(l-r)a_]2 (34) 
FTibv,bA,) = il-2ruj + r^)[iuj^-l)bl + bX] 
FL{bA„bA2,bA^) = [{uJ-r)bA, + iuj^ -l)irbA2+bA,)]^■ 
We compare our model results for the full rates Th with those of other models in (|35|), in 
units of 10^^|V^fep s~^. (Note that I+/2 and 1^^^ are linear combinations of ^Pi and ^Pi. 
Results computed in the latter basis are denoted by asterisks in (^), and "HQ" denotes 
heavy-quark limit results.) 
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N Jf 



3l 



^1/2 
-'-1/2 

0^/2 

1 + 
-■-1/2 

1 + 
-^3/2 

9+ 

^3/2 
I3/2 

^3/2 
0v2 

^1/2 



ISGW 
(full) 

[19] 



CNP 
(full) 



ours 
(full) 



ours 
(HQ) 



CNP 
(HQ) 

m 



1.25 

3.24 

0.049 

0.180* 

0.045* 

0.061 



0.013 
0.010 



0.76 

2.3 

0.076 

0.076* 



0.900 
2.91 
0.028 
0.040 
0.117 
0.076 
0.0021 
0.0003 
0.055 



0.771 

2.32 

0.013 

0.026 

0.070 

0.104 

0.0002 

0.0002 



0.723 

2.22 

0.026 

0.036 

0.052 

0.103 



Suzuki 
(HQ) 



0.741 

2.20 

0.011 

0.015 

0.034 

0.049 

0.0007 

0.0007 

0.045 

0.115 



(35) 



There are clear similarities among the various models. At the same time, measurements of 
these exclusive resonant processes could potentially discriminate between the models. 

In Fig. 6 we show the electron energy spectra for the various resonant final D states and 
their total contribution to the inclusive B decay spectrum. The location of the maximum 
agrees with ISGW . Fig. 7 shows the spectra of the excited D final states in more detail. 



Part III 

Some details of the heavy-quark limit. 

We first recover one of the basic properties of heavy quark effective theory. The following is 
satisfied identically: 

{P\Qix)it^ - mQ)Q{x)\P) = (36) 

when the pseudoscalar meson P is on shell. This is easily seen because — mq just removes 
one of the heavy-quark propagators from the graph for (p6D, leaving the P mass function 
evaluated at = M^. This vanishes, by the definition of the mass. In the heavy-quark 
effective theory at lowest order we have Q{x) ~ exp(— irngf ■ x)hy{x). Using '^h^ = h^, we 
find 

{P\Q{x){t^ - mQ)Q{x)\P) ~ {P\K{x)td ■ vK{x)\P) = 0. (37) 

Thus we recover the equation of motion iv ■ Dh^ = of heavy quark effective theory; here 
D = d since there are no gluons in the model. 

The heavy-quark limit results in ( pSj) are included for comparison purposes only. They 
have limited physical meaning since they are computed using the above formulae with Isgur- 
Wise functions for the form factors, while the meson masses are kept at their physical values. 

It has become standard to display the differential oo spectrum by taking the square root 
of (|32|) and dividing out model-independent factors in such a way as to make the result equal 
to the corresponding Isgur-Wise function if the form factors are replaced by their values in 
the heavy-quark limit. This procedure is applied to the 1/2+ and 3/2+ final states (the 
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curves labelled "full") in Figs. 8 and 9, respectively. We note the considerable deviations 
from the heavy-quark limit. 

An important quantity is the lowest-order mass difference A between meson and heavy 
quark: 

M = mQ + A + 0(AVmQ). (38) 
We extract the following approximate values for the = 1 states: 

j[ 1/2- 1/2+ 3/2+ 3/2- 



A 500 860 990 1330 



(39) 



These numbers make clear the breakdown of an expansion in A/rric for the higher-lying 
states. 

We also display the following zero-recoil values of the Isgur-Wise functions, along with 
our previous result for the slope of the usual Isgur-Wise function ^ = ^1/2 at zero recoil. 

n/2(l) r3/2(l) ^3/2(1) . . 

1.28 0.21 0.29 0.013 ^ ' 



In comparison, the ISGW model result is ti/2(1) = ^3/2(1) = 0.315 [g3|. The QCD sum rule 



approach of |^ found ti/2(1) — ^3/2(1) = 0.2 to 0.25]^ while the Bethe-Salpeter approach 



of [H found ri/2(l) = 0.21 and r3/2(l) = 0.44. 

Another issue is the extent to which the Bjorken sum rule is saturated by the states we 
have calculated. This sum rule reads 

- el/2(l) = 1/4 + |n/2(l)r + 2|r3/2(l)r + • • ■ (41) 

where the ellipsis denotes contributions from higher radial excitations, from states with 
quantum numbers other than 1/2", 1/2+ and 3/2+, and from inelastic continua. We find 
numerically 1.28 = 0.25 + 0.04 + 0.17 + . . . = 0.46 + ... for Eq. (|T]). The sum rule is far 
from being saturated by the resonances we have considered so far; this result reflects the 
situation already observed in the full model. As already stated, our model indicates that this 
is due mainly to nonresonant contributions to final states containing a pion. The discrepancy 
is even larger here because the sum rule was derived in the heavy-quark limit, and the D 
mesons are rather far from this limit in our model. 

One last issue of current interest is the dependence of the slope of the Isgur-Wise function 
on the light-quark mass. We find that the slope increases when the light-quark mass is 
replaced by the strange-quark mass, from -1.28 to -1.55. This behavior agrees with other 
analyses in quark models [^, sum rules [^] and lattice ^7j, but it disagrees with heavy- 
quark chiral perturbation theory [^. Our dependence of the slope on the light-quark mass 



is given in and it is far from linear. 



The decay D KHv 

We now explore how the model may be extended to describe strange mesons. We will 
consider the K* and will avoid the K and its complicating pseudo-Goldstone nature. The 

^The range of values reflects the conflict between Fig. 1 and Table 1 of 
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process K* (892)'^ e~^iye is of interest since models tend to disagree with the data by 

predicting a branching fraction which is larger than the measured value. For us, one option 
is to model the K* exactly like we model the D* and B* mesons. This would be making 
the rather suspect assumption that the s quark acts like a heavy quark. The other option, 
perhaps equally suspect, is to treat the s quark as light, on the same footing as the u and d 
quarks. This would entail modifying the vertex damping factor in in order to treat the 
d and s symmetrically. The damping factor then becomes 



-B + A2 

a 



+ 



+ A2 



(42) 



We will give the results for these two options. 

With a strange quark mass of 250 + 170 = 420 MeV, we find A^* = 598 MeV and 
Zk* = 836 MeV. The relevant observables can be compared to the measured values as 



follows. 10 



T 



r+/r_ B{D 



+ 



K 



*o 



experiment 
model (heavy s) 
model (light s) 



1.23 ±0.13 
0.99 
1.04 



0.16 ±0.04 
0.24 
0.25 



4.8 ±0.4 
8.1 
2.8 



(43) 



We see that only the branching fraction is very sensitive to how the s quark is treated, and 
interestingly enough, the experimental branching fraction lies between the two extreme ways 
of treating the s. 

We may also extract the form factors V{q'^), Ailq"^) and ^2(1?^) for q'^=0, which may then 
be compared to nonrelativistic quark model and to sum rule results. For additional models 

ig, 

^(0) Ai(0) ^2(0) 



see 



model (heavy s) 
model (light s) 
ISGW 
BBD 
experiment 



1.14 
0.68 
1.1 
1.1 
1.1 ±0.2 



0.76 
0.48 

0.8 

0.5 
0.56 ±0.04 



0.82 
0.52 

0.8 

0.6 
0.40 ± 0.08 



(44) 



Note that the experimental values here |TD| have been extracted by assuming specific pole 
expressions for the form factors. We may compare those form factors to our model form 
factors by fitting our form factors in the physical region to the pole form F{q^) = F(0)/(1 - 
q^ /my). The values for mi? in GeV are the following. 





V 


Ai 


A2 


assumed in [10] 


2.1 


2.5 


2.5 


model (heavy s) 


1.65 


2.3 


1.8 


model (light s) 


2.0 


3.8 


2.4 



(45) 
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The form of the vertex damping factor 



The model relies on a vertex damping factor to suppress the flow of large Euclidean momenta 
through light-quark propagators, and a specific pole form has been chosen. 

Z2 



Fie 



(46) 



-fc2 + A2 

This introduces singularities into our quark loop amplitudes for momenta outside the physical 
region P] . Viewing our amplitudes as functions of the quark masses while keeping the meson 
masses fixed, there are also cusp-like singularities for special (presumably unphysical) values 
of the quark masses 0. We wish to explore here how our conclusions depend on our choice 
of the damping factor, and on the singularity structure of the model in particular. To do 
this we will substitute an exponential form for the damping factor, 

A2 



Rexp { — 



(47) 



This may be less realistic for large — fc^, but it is free from unphysical singularities at finite 
— fc2. The parameters R and A will be determined for each meson from the self-energy graphs, 
as was done for the quantities Z and A of (0). Our standard quark masses imply Rb=3.81, 
Ab=Q20 MeV, i?D=2.37, Ad=518 MeV, Rd*=2M and A^.=785 MeV. The shape of the 
exponential damping factor is compared to the original damping factor for the B meson in 
Fig. 10. 

We give the values of the two form factors for B D^*Hv which at u!=l are protected 
from first order corrections in the heavy quark expansion. 



vertex 


rric 


bAj2 


a+ 


original 


1.44 


1.16 


1.11 




1.45 


1.12 


1.17 


exponential 


1.55 


.95 


1.83 




1.50 


.97 


1.34 




1.44 


1.01 


1.19 




1.40 


1.06 


1.17 




1.35 


1.15 


1.17 




1.30 


1.56 


1.19 



(4J 



As we have explained elsewhere 0, the unphysical procedure of changing the quark masses 
while holding the meson masses fixed can lead to clearly unphysical values for amplitudes. 
The interesting result in the exponential case is the very large corrections to a+ for any 
reasonable value of rric. Thus our conclusion regarding the existence of large corrections to 
the heavy-quark limit in at least one of i? ^ D*£v oi B ^ DiV H remains unchanged. 



Conclusions 

We have thoroughly explored a relativistic quark model for hadrons containing one heavy 
quark. It is extremely economical in its description of nonperturbative QCD, its only param- 
eters being the quark masses, yet it is consistent with the symmetries of QCD. It displays 
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overall agreement, across a broad range of different processes, with the data when it exists 
and with more complicated models such as QCD sum rules in the absence of data. The main 
difference with other models is the suggestion that large nonperturbative departures from 
the heavy quark limit occur in certain quantities, most notably in the processes 
and in the heavy-quark contribution to the DD*^ coupling. At the same time, it is consis- 
tent with data for the shape of the spectrum for B D*lv and for the branching ratios for 
D* —>■ D'-f. To test the pattern of deviations from the heavy-quark limit in our model and 
others, we especially encourage efforts to extract and compare Vet from each of -B ^ D*lv, 
B Diu and ^ AJiy. 
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Figure captions 



FIG. 1: Comparison of model prediction for B 
notation is the same as in M. 



D*iV with recent CLEO data |[. The 



and flSOD. Certain coefficients 



FIG. 2: Full-model form factors q and di, defined in 
have been divided out as explained in the text. Also shown is the Isgur-Wise function ri/2, 
defined in (0). 



FIG. 3: Full-model form factors fi and (^j, defined in (^Uj) and (PPD- Certain coefficients 
have been divided out as explained in the text. Also shown is the Isgur-Wise function t^/2i 



defined in 



FIG. 4: Full-model form factors hi and fcj, defined in (|29|) and (|30|) . Certain coefficients 
have been divided out as explained in the text. Also shown is the Isgur-Wise function ,^3/2, 
defined in 



FIG. 5: Full- model form factors for B Diu. a± are for the ground-state D, and l± are 



for its first radial excitation, as defined in (|29D and (|30|). 



FIG. 6: Electron energy spectra for resonant semileptonic B decays in the full model. 

FIG. 7: Electron energy spectra for resonant semileptonic B decays to excited D states 
only. 

FIG. 8: Differential spectrum for i? ^ 1 l^^iu. Normalization is such that the spectrum 
becomes equal to the Isgur-Wise function T1/2 when the form factors are replaced by their 
values in the heavy-quark limit. 

FIG. 9: Differential spectrum for B ^ 1 Xt^i^iv- Normalization is such that the spectrum 
becomes equal to the Isgur-Wise function r3/2 when the form factors are replaced by their 
values in the heavy-quark limit. 

FIG. 10: Comparison of usual damping factor with exponential damping factor for B meson, 
with physical values of parameters. 
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